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The ATLAS experiment has recently observed a 3cr excess in a channel with a leptonically decaying 
Z, jets, and J?“ lss . It is tantalizing to interpret the signal as the first sign of a natural supersymmetric 
spectrum. We study such a possibility in a minimal model containing light stops and a neutralino 
LSP. The signal is characterized by a novel topology (compared to previous attempts) where the 
Z is emitted from a colored particle in the first step of a decay chain, namely f 2 —► tiZ which 
is characteristic of mixed stops. We show that the excess is compatible with a compressed stop 
spectrum and is not excluded by any other relevant search, finding some regions of parameter space 
with signal strength within lcr of that measured by the ATLAS collaboration. In addition, we notice 
that the corresponding CMS search could be prone to background contamination in unexpected 
topologies of this kind. 


I. INTRODUCTION 

Supersymmetry (SUSY) is a leading candidate for 
resolving the large hierarchy problem of the Stan¬ 
dard Model. In the Minimal Supersymmetric Standard 
Model (MSSM) and simple extensions, a necessary fea¬ 
ture for a complete resolution of the hierarchy problem 
is the presence of two light (sub-TeV) colored stops and 
one light left handed sbottom (to accompany the left 
handed stop). A common assumption in these models is 
an exact R-Parity, and the presence of a neutral, stable 
lightest Supersymmetric particle (LSP). In this case, if 
the third generation squarks are accompanied by a neu¬ 
tralino LSP, x°, then typical decays of these particles 
include U ,2 —> t^x° (where the superscript on t indi¬ 
cates the possibility that it is off-shell), t 2 —t I\Z and 
bi —> bx°. The signatures of this scenario are therefore 
jets, missing transverse momentum (U“ lss ), leptons and 
^-tagged jets. 

Dedicated searches for 3rd generation squarks have 
found no deviations from SM predictions, placing strin¬ 
gent constraints on its parameter space. On the other 
hand there remain significant windows allowing the 
mass of the lightest stop toj to be as light as 200 GeV, 
provided that there is a compressed spectrum which 
softens the px distributions of the final state particles. 
Intriguingly, a recent ATLAS search found a 3cr excess 
in final states containing a leptonically decaying Z bo¬ 
son, jets, and large £vj? lss [Tj. They found 29 events in a 
combined signal region with expected SM background 
of 10.6 ± 3.2 events. We wish to explore the possibility 
that this excess is a first signal for direct production of 
t 2 followed by the decay t 2 —>• U €1 
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1 See |2] for other recent work on this signature 


Various attempts have been made to explain this 
excess in terms of SUSY modelij^] R1TH?I. In all of 
these studies, pair produced colored particles (squarks 
or gluinos) decay into quarks and an uncolored particle, 
which then decays into a Z boson and an LSP. The prin¬ 
cipal challenge faced by these models is in explaining the 
ATLAS excess while simultaneously evading the many 
bounds imposed by other searches by ATLAS and CMS 
for multileptons or jets and U™ lss , as well as a similar 
CMS search for the same final state m which saw no 
excess over Standard Model (SM) backgrounds. This 
latter search imposed different cuts from the ATLAS 
search and so does not necessarily rule out new physics 
explanations for the ATLAS excess, yet it still imposes 
stringent constraints (see, e.g. !)■ 

The phenomenology of the signal proposed in this pa¬ 
per differs from that of the aforementioned possibilities 
in several key respects. Firstly, the topology differs in 
that the Z boson is emitted at the first stage in the 
decay, rather than at the end with the LSP. This opens 
up the possibility that the CMS search is subject to 
significant background contamination, as we discuss at 
the end of section m More significantly, our scenario 
requires the presence of three new colored particles in 
the spectrum which are lighter than in previous expla¬ 
nations, the heaviest of which gives rise to the desired 
signature. Evading dedicated searches for these parti¬ 
cles places very particular constraints on the mass split¬ 
tings and decays of the squarks. As we shall discuss in 
section [TT1 this requires a compressed splitting between 
ti and x> an d possibly also between t\ and t 2 . This 
in turn motivates the consideration of flavor violating 
decays of t\ into ux° or cy°, resulting from mixing be¬ 
tween the right handed squarks. Such mixings have 
been discussed in recent years motivated by the ques- 


2 See [3] for a discussion of this excess in the framework of Com¬ 
posite Higgs models. 




tion of natural SUSY and light stops HMS], but with¬ 
out the Z decay necessary to explain this excess. 

In section[TT]we provide a systematic discussion of the 
possibilities that this minimal stop scenario affords for 
explaining the excess, identifying three distinct scenar¬ 
ios characterized by the mass splittings involved and 
the assumed decay mode for the light stop. We pro¬ 
ceed in section |III| to describe the main experimental 
searches placing limits on these scenarios, and perform 
scans of their parameter spaces to find regions in which 
the excess can be explained while evading those limits in 
section |TV} We also note the possibility that the signal 
topologies that we have identified could cause significant 
background contamination in the CMS on -Z search. 

II. MODEL OVERVIEW 

In this paper we assume a minimal model including 
two light stops, <1 and and a left-handed sbottom 
b\. We also require one neutralino LSP, x°. As will 
be discussed below, the identity of this neutralino is 
not relevant to collider phenomenology for two of the 
scenarios that we consider, while for the third case it 
will be assumed to be mostly Binc0 All other SUSY 
states are assumed to be heavier than these particles 
and decoupled. In the MSSM, the stop and sbottom 
mass matrices in the gauge-eigenbasis are given by m 


basis 



The sbottom mass is then given by 

= m h c l- + TO I a 0 { - m t~ A ul + A d L ■ ( IL4 ) 

We assume the decoupling limit for the Higgs sector, 
so that the Higgs mixing angle a = /3 — 7r/2. The phe¬ 
nomenology of this simplified model varies only slightly 
with tan (3 and we therefore choose to fix tan/3 = 20. 
The remaining free parameters in the model are , 
nif , o, and cg i . Even with such a modest amount 
of new particles, this model admits a rich phenomenol¬ 
ogy with many possible final states, depending mainly 
on the assumed mass splittings and mixings involved. 
We seek scenarios with a large branching ratio (BR) for 
t 2 —> t\Z, and which are poorly constrained by dedi¬ 
cated searches for t-\ and b\. In the following subsec¬ 
tions, we systematically discuss the various possibili¬ 
ties and present a categorization of interesting scenarios 
based on the assumptions made about the decays of t\ 
and the mass splitting between t\ and t 2 - 

A. ti decays 


n? = ( + m t+ v ( a t s 0 - mcp) \ 
* \v (a t sp - H*ytCp) ml 3 + ml + A u R ) 

m ? = f m 2 Q3 + A Jl v (a* b c 0 - yy b sp)\ 
b \y ( a b C 0 - y*y b sp) m 2 d3 + A Sr J 


(H.1) 

(H.2) 


where itiq 3 , m U3 , mj 3 , a t , y are soft SUSY breaking 
parameters, cp and sp denote the cosine and sine of /3, 
and Aq = (T^ — Qq sin 2 9 w)c 2 pm 2 z with T$q and Qq 
denoting the third component of weak isospin and elec¬ 
tric charge respectively. The Higgs vev v is ss 174GeV. 
We assume the right handed sbottom is decoupled, with 
m 2 - rri'n . We replace the MSSM parameters in these 

CL 3 ' £ e 3 

mass matrices with physical parameters: the stop mass 
eigenstates m bi , , and the mixing angle 0 < Of < 7r/2 

which rotates the gauge-eigenstate basis into the mass 


! We do not want to address any cosmological issue in this work 
but let us notice that a stable Bino is overproduced in the early 
universe according to the usual thermal freeze-out calculation. 
On the other hand, a small ti— x° mass splitting of 0(30 GeV 
allows for the possibility that the correct relic density results 
from stop-neutralino coannihilation 11711181 . Alternative solu¬ 
tions simple solution are to assume a low reheating temperature 
or that that the Bino is actually the NLSP (for instance with 
a gravitino LSP) but still stable on detector lengths. 


The strongest constraints on the t,-\ apply if it decays 
directly to a neutralino and on-shell top, leading to final 
states with large E^ lss , hard &-jets, and leptons. We 
therefore take the splitting m bl — m^o < mt, such that 
the only flavor-conserving decays that are kinematically 
available to the light stop are into the three- or four- 
body final states Wb\° or ff'bx 0 (where //' are pairs 
of fermions that may be produced in the decay of an off- 
shell W). This allows t\ to be as light as 300 GeV for 
generic values of this splitting, and as low as 200 GeV 
in some narrow windows of parameter space (see [20] 
for a detailed discussion). 

Due to the substantial kinematic suppression of the 
partial width into these states, it is possible that fla¬ 
vor violating decays might dominate even with small 
couplings. This motivates our consideration of flavor 
violating decays. One well explored possibility arises 
even with Minimal Flavor Violation (MFV) [2T|, in 
which case it is possible that loop-induced decays into 
charm and neutralino can dominate over four-body de¬ 
cays (22]. In recent years an alternative scenario has 
been explored, that non-MFV mixings between right 
handed up-type squarks can substantially alter stop 
phenomenology. Briefly, the essential point for our 
analysis is that the strongest constraints on the size of 
squark flavor mixings from low energy observables ap¬ 
ply to the down sector, and on mixings between up and 
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FIG. 1: Flavor conserving (left) and flavor violating (right) decays contributing to the on -Z excess. In the flavor 

violating case, q can be either an up or a charm quark. 


charm squarks. The constraints on the down sector also 
impose constraints on the left-handed up type squarks. 
Crucially, there are no direct constraints on £r — cr or 
tR — ur mixings individually, but only on their product 
(coming from the D° — D° system). We refer readers to 
the papers mm for a more detailed discussion. As a 
consequence, there may be size-able mixing between tR 
and cr or tR and ur, but not both. 

The degree of flavor mixing can be parameterized by 
the quantity e = (tti?)* 3 / (m\)u, where m~ is the up- 
type squark mass matrix in the Super-CKM basis and 
i is 1 or 2. We do not require 0(1) mixings in order 
to change the decay patterns of the lightest stop, so 
long as it has an 0(1) admixture of tR. In particular, 
in the four-body regiorj^J e > 10~ 3 is sufficient for the 
decay i-\ —>■ qx° to occur at least 90% of the time, while 
e > 10~ 2 is sufficient for much of the three-body region. 
These flavor-mixing angles are sufficiently small not to 
play a noticeable role in the phenomenology of the i 2 , 
and are relevant for t\ only because its flavor-conserving 
decays are heavily suppressed. It can also be assumed 
that (m|)jj (m|) 33 , such that despite introducing 

a small admixture of first or second generation squark 
flavor into the two dominantly stop mass eigenstates, 
the other mass eigenstates are beyond reach of the first 
run of the LHC. 

We therefore consider separately scenarios where ti 
undergoes a flavor-conserving decay (denoted Fq ), or 
a flavor violating one (^). In either case we will assume 
a 100% BR for the light stop for simplicity, though as 


4 Even when talking about flavor violating two-body decays, we 
find it convenient to label the regions of parameter space in 
the ii, x° plane by the possible flavor conserving decays. The 
‘three-body region’ is defined by myy - l-mj, < — m~o < mt, 

while the ‘four-body region’ is defined by mj < m 7 | — 771 -0 < 
m w + m b . 


pointed out in ;23j the limits on light stops can be sub¬ 
stantially reduced for mixed ^ and Fq decays. The fla¬ 
vor violating decay may be into u\° or cry 0 , but not 
both. The only difference this makes regarding collider 
phenomenology is that the constraints on final states 
containing charm quarks are often stronger than on up 
quarks, due to significant progress made on charm fla¬ 
vor tagging by the experimental collaborations. This 
will be discussed in more detail in section m It should 
be noted that the precise measurement of the neutron 
Electric Dipole Moment (EDM) places constraints on 
tR ~ ur mixing in the presence of large stop L-R mixing 
as exists in our model. This comes from loop contribu¬ 
tions to the up quark EDM involving gluinos, and de¬ 
pends sensitively on details about the particles which we 
have assumed to be decoupled in our scenario. Nonethe¬ 
less, it is demonstrated in [24] that e < 10 -2 can be 
consistent with the EDM constraints without making 
unnatural assumptions about the masses of the other 
particles, or about cancellations between different con¬ 
tributions. 


B. t2 decays 

The second important distinction to be made between 
different classes of scenarios relates to the mass split¬ 
ting between t 2 and t,\, and the role that this plays in 
determining the branching ratios of the three squarks. 
For sufficiently large mass splittings, the possible two- 
body decays of the heavy stop are £ 2 -» £1 Z, t 2 —> t-\ h, 
£2 —> biW + and £ 2 —► tx°. The BRs into these states 
are most sensitive to the mixing angle cg r and for split¬ 
ting 150 GeV < rnf 2 — < 300 GeV the BR into Zt\ 

is maximized at a value between 0.6 < BR(f 2 —t t.±Z) < 
0.8 for 0.5 < cg { < 0.55, as illustrated in Fig. ([ 2 ]). Since 
we are interested in maximizing the f 2 —> t\Z BR we 
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FIG. 2: Heavy stop branching ratios in the split 
scenario, for mg 2 = 450 GeV, = 250 GeV, 
myo = 210 GeV, tan/3 = 20. 


can use this to fix the mixing angle. This class of sce¬ 
nario, which we label ‘split’, has a three-dimensional 
parameter space in m g , mg 2 and myo. Note that for 
small cos 9 t , t\ and b± are almost degenerate and hence 
the decay t\ —> b\W is not kinematically available. 

An interesting alternative is that the £2 — ti mass 
splitting is sufficiently small that the only two-body 
decay kinematically allowed for t 2 is t 2 t\Z. For 
—mi 2 < rrih, the decay t 2 —> t\h is forbidden. For a 
wide range of cg { , the decay t 2 —► b±W + is not kinemat¬ 
ically available. Finally, if +7712 < mg 2 < m t +my 0 , 
the only two-body decay that is kinematically available 
is the Z decay. Combining the bounds gives the re¬ 
quirement that — myo < mt — mz ~ 85 GeV. 
Coincidentally, this turns out to overlap almost exclu¬ 
sively with the four-body region, which is defined by 
— myo < mw + nib ~ 85 GeV. This condition de¬ 
fines what we call the ‘compressed’ scenarios, in which 
the only two-body decay available to the heavy stop is 
into t\Z. The t 2 BRs are therefore insensitive to cg- 
in this regime. This parameter does control the mass, 
and therefore also the decays of 61 . A heavy bi can de¬ 
cay into t-iW~, but if this is not available it will have 
~ 100% BR into fey 0 . For generic mass splittings this 
decay is highly constrained by dedicated searches, but 
having chosen small splittings for t\ — x° and t 2 — t\, 
this automatically places also the sbottom decay into 
the compressed regime in which this channel is particu¬ 
larly challenging for experimental searches. As we shall 
discuss in more detail in section m for a broad range 
of stop mixing angles the sbottom lies in a funnel of 
parameter space not constrained by these searches. 


Scenario 

ti decay 

m to - m u 

BR(*2 —^ tiZ') 

F c -C 

ffbx° 

100 GeV 

1 

f-C 

cx° / ux° 

100 GeV 

1 

p-s 

cx° / ux° (125 - 300) GeV 

0.7 ± 0.1 


TABLE I: The three scenarios considered in this 
paper, labeled Flavor Conserving Compressed, Flavor 
Violating Compressed, and Flavor Violating Split. 
The Fq /f designation refers to the decays of t-\ , and 
the compressed/split designation refers to the splitting 
between t 2 and t\. 


C. Three Scenarios 

Four combinations of Fq versus f and split versus 
compressed are possible. Flavor conserving split (Fq - 
S) is the most highly constrained by dedicated searches, 
and we have been unable to find a region of its parame¬ 
ter space which permits an explanation of the Z excess 
without being excluded by other searches. We there¬ 
fore do not discuss this possibility in detail in this pa¬ 
per. Three combinations remain, which are summarized 
in table (jl|). For the compressed scenarios, we choose 
mg 2 — =100 GeV and explore the — myo plane. 

We set the branching ratios BR(i 2 —> h X°) = 1 and 
BR(&i —> bx°) = 1. In principle, t 2 also has compet¬ 
ing three-body decays to b±, but these are sensitive to 
mg and would have a BR no more than 0(5%). We 
neglect this effect for simplicity. For the split scenario, 
we explore t 2 — 1\ mass splittings between 125 and 300 
GeV, choosing cg f = 0.5 which is close to the optimal 
value for maximizing BR(t 2 —> t\Z) over most of the 
parameter space. This angle then also sets the BRs for 
the b-\ . We also compute the t 2 -7- tx° and b\ —> bx° 
BRs assuming a Bino LSP. 


III. RELEVANT SEARCHES 

ATLAS and CMS have a wealth of searches looking 
for large MET with all types of additional particles in 
the final state, each potentially providing a limit on stop 
and sbottom production. Since in general ti,t 2 , and b 1 
will contribute to each bound, the constraints needs to 
be recast with care. Our goal is to examine all the pa¬ 
rameter space with the simplified topology discussed in 
section [IT] For the compressed cases this involves a scan 
in the m^o — plane while for the split case it in¬ 
volves a scan in both my 0 — mg and mg x — mg 2 . We will 
be exploring scenarios with compressed mass splittings, 
especially between t\ and \ 0 . This is a region that 
is very challenging experimentally. The most robust 
and model-independent limits in the most compressed 
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regime come from dedicated searches for events with 
a hard jet coming from Initial State Radiation (ISR), 
which do not depend sensitively on the details of the 
£-1 decay. Searches for jets and U™ lss are highly con¬ 
straining for spectra with a large splitting between the 
LSP and colored particles, but are challenging to in¬ 
terpret in the compressed regime where there are large 
systematic uncertainties. Searches involving 6-tagged 
jets place important limits on some of our decay chan¬ 
nels. Finally, there are dedicated searches for events 
containing Z bosons which could be sensitive to our 
model. We discuss the details of these searches in the 
following subsections, beginning with the ATLAS on-Z 
search with a 3<r excess. 


A. ATLAS on -Z 

The ATLAS on -Z search looked for final states with 
two leptons with invariant mass around the Z- pole, 
Ef ss > 225 GeV, and H T = + + Vt + Eiejete Pt > 
600 GeV. The Ht and E™ lss cuts pick out events 
with large kinetic energies. ATLAS found 16 events 
(4.2 ± 1.6 expected) in the electron channel and 13 
events (6.4±2.2 expected) in the muon channel for a to¬ 
tal of 29 events (10.6 ± 3.2 expected). Running pseudo¬ 
experiments they concluded that this corresponds to a 
3.0cr deviations from the SM. 

To estimate the number of events needed to explain 
the signal we use a log-likelihood method and profile 
over the background uncertainties using a Gaussian ap¬ 
proximation. Using asyptotic formulae (25] to estab¬ 
lish two-sided convidence limits (CL), we find that a 
minimum of 7.1 (12.4) signal events are required to be 
consistent with the excess at the 95% (68%) CL. 


scribed in section IV) We found that the mid E™ lss 
bin usually provides the dominant limit, but the high 
E™ lss bin is sometimes competitive. The low E™ lss bin 
is never competitive with the other two in our simula¬ 
tions 0 


c. —y tiZ 

Searches for the signal t 2 —=► t\Z at 8 TeV have been 
performed both by ATLAS (21] and CMS .28]. The 
searches have competitive bounds, with ATLAS hav¬ 
ing a slightly better exclusion. For simplicity we only 
use the ATLAS search to place bounds on the Fc- 
C scenario. The main cuts in this search are on the 
invariant mass of the leptons (which are required to be 
around the Z pole), the number of jets, number of lep¬ 
tons, and the pt of the reconstructed Z. We find that 
the reach of the search is limited in the compressed 
regime, both due to the small t 2 ~ti splitting and the 
small ti~x° splitting. Firstly, the search regions requir¬ 
ing two leptons require a boosted Z candidate which 
is suppressed in the Fq-C scenario by the small G+ 
splitting. Secondly, we find that the soft leptons and 
6-jets coming from the decay of t\ in the 4-body regime 
reduce the acceptance in the 3-lepton bins and the 6- 
tagging efficiency. Having recasted this search for the 
Fq-C scenario, we find that it does not place competi¬ 
tive limits and we therefore omit it from our scans. We 
note that the CMS search is optimised for mass spitting 
—m ^o ~ nit and has potential sensitivity to our Fc- 
C scenario only in its 3-lepton bins, and will therefore 
also have degraded sensitivity in the 4-body region. 


D. Jets+MET+0/1 lepton 


B. CMS on -Z 

CMS performed a search analogous to that done by 
ATLAS, looking for events with opposite-sign same- 
flavor leptons and E™ lss PS] which provides an impor¬ 
tant constraint on our model. CMS split their on -Z sig¬ 
nal regions into six bins, depending on jet multiplicity 
and E™ lss . Three bins measure events with nj e ts > 2, 
and three use nj ets > 3. Our signal model does not 
produce a significant number of 2-jet events, and we 
therefore place constraints on our scenarios using the 
3-jet inclusive bins which have higher expected sensi¬ 
tivity. The three bins in this category are split into the 
following Ef iss windows: 100 GeV < Ef iss < 200 GeV, 
200 GeV < £+ ss < 300 GeV, Ef iss > 300 GeV. We 
use the most constraining of these bins to constrain our 
scenarios at each point in parameter space, using the 
95% CL S limit to set bounds [25]. In our scans (de- 


The MSSM-inspired jets+MET searched provide an 
important constraint on our scenarios as ti,t 2 , and 
&! can all contribute to this signal. ATLAS 22] and 
CMS (20: have both performed searches for this sig¬ 
nature at 8 TeV and interpreted them in terms of a 
variety of SUSY models, including direct production of 
squarks decaying via q —> jx° (which is identical to 


5 A possible concern is that the combined limit from different 
bins could be more severely constraining. We find that usually 
only one bin is constraining and we do not expect a combination 
to notably alter the limits. 

6 It was pointed out in m that searches using shape-based anal¬ 
yses might have better sensitivity for ti —> jx° then jets+MET 
in the limit of small — m^o. The bounds were computed 
using 7 TeV data and are not constraining compared to the 
8 TeV jets+MET search. It would be interesting to see how 
these would change with the full 8 TeV data set. 
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t-\ in the f scenarios). In this work we focus on the 
200 GeV < < 350 GeV region but CMS only 

presents limits for to, > 300 GeV. Thus to study the 
bounds from these searches we use the ATLAS analysis. 
Jets + MET searches in these regions are particularly 
challenging due to the low pp of the outgoing parti¬ 
cles. As a consequence their search does not constrain 
single squark production decaying in this channel for 
rriy-o > 160 GeV. Another key factor which limits the 
reach of this search in the compressed regime is the 
systematic uncertainty on the (acceptance x efficiency) 
associated with the high sensitivity to ISR. ATLAS pro¬ 
vides uncertainties for each signal region in their aux¬ 
iliary material (available on the ATLAS public results 
website), and these range from 10% to 50% in constrain¬ 
ing bins. 

In order to interpret the results of this search in terms 
of limits on our scenarios, it is necessary to combine 
the contributions from t 2 , t\ and b±. Each of these 
channels will come with its own set of uncertainties 
which vary from bin-to-bin. Without a dedicated de¬ 
tector simulation it is not possible to robustly account 
for these effects. Instead, we estimate reasonable sizes 
for these uncertainties and study the effects of vary¬ 
ing these assumptions. The ATLAS collaboration also 
interpreted their results in terms of some multi-step de¬ 
cay chains, and find uncertainties that range from 10% 
to 80% in similarly compressed regimes. Using the ‘r* 
methocQ 3T|, we find good agreement with the ATLAS 
exclusion on t\ if we assign a uniform systematic un¬ 
certainty of 30% on its signal strength across all signal 
regions We then assign a nominal 30% uncertainty also 
on the t 2 and 6j contributions, consistent with the afore¬ 
mentioned uncertainties quoted by the ATLAS collabo¬ 
ration for other compressed multi-step decay processes. 
In order to asses the sensitivity of our results to this 
choice, we also vary the uncertainty on the t 2 and b\ 
channels to 20% and 40%. We used ClieckMATE [5T 
to apply this technique to all our scenarios however we 
found the bounds of this search to be subdominant in 
all cases except for the JP-S case. This is because this 
is the only scenario in which we venture close to the 
existing bounds on q —> jx°, where a combination with 
the other channels might then result in an exclusion. 

In addition to jets+MET there are searches which re¬ 
quire an additional isolated lepton by both ATLAS [32] 
and CMS [33], which are potentially sensitive to the 
Fc-C scenario. However, the limits on the light stop 
are weaker than the other limits which we consider for 


7 In the r-method, a signal is excluded if r = (S 
1.96AS)/Sgg s > 1, where AS is the systematic uncertainty on 
the signal strength and Sgg s is the limit on the signal strength 
at the 95% confidence limit. 


this search, and we find that the heavy stop produc¬ 
tion does not contribute significantly to this signal. We 
therefore do not include this limit in our scans. 


E. Sbottom bounds 


The sbottom mass is determined by the mixing angle 
as given in eq. (II.4). The phenomenology of the sbot¬ 
tom differs substantially between the compressed and 
split scenarios. In the compressed scenario, the only 
two-body final state available for it to decay into is 6y°. 
Furthermore, the mixing angle is a free parameter since 
the branching ratio of t 2 —> t\Z is fixed to 1 by the 
kinematics. On the other hand, in the split scenario 
this mixing angle is fixed at cg f = 0.5 and the sbot¬ 


tom decays almost exclusively into i\ W in most of the 
parameter space. 

The bi —» b\ 0 decay of the compressed scenar¬ 
ios is constrained by dedicated CMS and ATLAS 
searches [Ml (35], and we focus on the ATLAS analy¬ 
sis because it places stronger limits in the compressed 
regime. This search places strong limits on this chan¬ 
nel, but allows for a sbottom with mass mr ~ 250 GeV 
if it has a small mass splitting with the neutralino. We 
also require that the sbottom is heavy enough to for¬ 
bid the decay t 2 —> biW, i.e., my 2 < mj >i + mw or in 


terms of t 2 — t\ mass splitting, +A m^ — mw- 

Combining the bounds gives 


< m' im (?n^o) — Am t -+ mw (III.l) 

where is the maximum allowed value of the 

sbottom mass for each m^o, coming from the ATLAS 
limit. This provides an additional constraint on the 
possible values in the m^o , plane. The additional t \, 
t 2 production channels are not expected to contribute 
significantly to this search. 

In the split scenario, ^-S , the sbottom decays pre¬ 
dominantly in b\ —^t\W — > jx°W. While there are 
no direct searches for this signal, there are searches for 
q — > x^W -» jx°W. We have checked the constraints 
due to this signal and found that we are well within 
experimental bounds for all regions of parameter space. 


F. Single high pr jet + 0,1,2 lepton 


In this work we are primarily interested in regions 
of parameter space with small t\ — x° splitting where 
many jets may not pass the pr cuts. In this case 
we have additional constraints from monojet searches. 
This search has been done by ATLAS for both t\ —> jx° 
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and ti —> bffW stop decay modes in a search by AT¬ 
LAS [86] . For the split case t\ is the only production 
channel which can replicate this signal, but in the com¬ 
pressed cases the b\ —>■ b\° could provide additional 
monojet events. To this end we recast the search in¬ 
cluding just ti production and both t\ and b i but we 
found comparable exclusions. For this reason we simply 
include the constraints computed by ATLAS directly in 
our analysis. 

CMS has performed a search for events with a high 
momentum ISR jet with the additional requirement of 
one or two soft leptons The preliminary results of 
this search provide the strongest existing constraints on 
the 4-body region of flavor-conserving t\ decays, sen¬ 
sitive to stop masses up to 320 GeV. The strongest 
bounds are derived from the 2-lepton signal region, and 
we therefore expect this limit to be highly sensitive to 
the i\ BR. 


G. Charm-tagging 

A final constraint on our signals are charm-tagging 
searches. ATLAS has two searches that employ charm¬ 
tagging looking for both i-\ —> c%° |36 as well as c —> 
C X° ES ■ The fi — > cx° search assumes the stop is in the 
four-body regime and hence is optimized for our signal. 
For this reason it has better sensitivity in our region of 
interest. For this reason we omit it from our plots. 

The c-tagging searches put constraints on f models 
which involve charm. For the p-S split scenatrio we 
find that these constraints rule out the region preferred 
by the ATLAS on -Z excess at the 95% CL. For this 
reason in this scenario we assume t\ decays to ux°- 
For the ^-C scenario the constraints are milder since 
we are exploring relatively large values. Recast¬ 
ing the ti —> cx° search we find that the limits on this 
scenario comparable to those directly on t\ alone, and 
thus we use the limits on this channel that are provided 
in [ggj . 


IV. SCAN 

For the scan we use Madgraph 5 v2.2.3 [39j, Pythia 
6.4 ||0] , and PGS [JT], including 1-jet matching. For 
jet clustering we use anti-kT algorithm with A R = 0.4. 
To roughly account for next-to-leading order (NLO) ef¬ 
fects we rescale our cross sections to their NLO values 
calculated by the SUSY Cross Sections group @3]. For 
jets+MET and double-checking monojet constraints we 
use CheckMATE [31], which makes use of DELPHES 
3 |43j, Fast Jet [44], and the anti-kT clustering algo¬ 
rithm [35] . For the compressed scenarios, Fc-C and 
p-C , there is one less free parameter (since the range 


Benchmark 

"A, [GeV] 

m h [GeV] 

m^o [GeV] 

p-value 

Fc-C 

380 

280 

260 

0.095 

f-C 

370 

275 

255 

0.17 

p-s 

450 

250 

210 

0.055 


TABLE II: Benchmark points chosen from the three 
scenarios. All other parameters are as described in 
section [n] The two-tailed p -values are calculated as 
described in section III A and a p -value of 1 would 
represent perfect agreement with the measured total 
event rate. 


of W( 2 has a relatively small range of viable options). 
We can perform a two dimensional scan over — m^o. 
p-S requires a 3 dimensional scan but for simplicity we 
scan over two slices in the parameter space. 

The scan showing the signal as well as limits from the 
different searches is shown in Fig. |3|, with the regions 
preferred by the ATLAS on -Z excess indicated by green 
shading and the contours labelling the 90% and 68% 
two-sided confidence intervals. The constraining 95% 
confidence intervals discussed in section |m] are shown 
by solid lines. The dashed line in the JAC scan indicates 
the limit on the decay ti -+ cy°, though the alternate 
decay t-\ — > ux° is also possible and not constrained 
by this line. The CMS on -Z limit in the ^-S scenarios 
is dashed as it has been calculated assuming no back¬ 
ground contamination. As we shall discuss below, con¬ 
siderable background contamination is expected which 
severely limits the sensitivity of the CMS search to this 
scenario. The jets+MET limit in the J?-S scenario is 
plotted with a band indicating the large uncertainties 
associated with this search in the compressed regime, 
as discussed in section EH The central line assumes a 
systematic uncertainty on the signal from all channels 
and in all bins of 30%. The band is obtained by varying 
the uncertainty on the t\ and bi production channels to 
20% and 40%. 

These plots indicate that all three scenarios can be 
consistent with the ATLAS on -Z excess at the 90% level 
and the two compressed scenarios at the la level, allow¬ 
ing for as many as 14 signal events. From each scenario 
we have chosen a benchmark point indicated by a black 
star in Fig. 0 , and detailed in table 0 - 


A. Kinematic distributions 

To further check the consistency with the data we 
compare our mu, E™ lss , Ht , and njets distributions with 
those measured by the ATLAS on -Z search (these cor¬ 
respond to Fig.6 and 7 in ID)- To compare the qual¬ 
ity of our signal we reproduce these plots in Fig. 0 
using the results by ATLAS to retrieve the SM back- 
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FIG. 3: Scans of the three scenarios. The regions favoured by the ATLAS on -Z excess are shaded green, with 
contours indicating 95% and 68% confidence intervals. Additional solid lines indicate the 95% limits described in 
the text. The dashed lines indicate limits under specific model assumptions that do not necessarily apply, as 
described in the text. The band on the jets+MET limit illustrates the considerable uncertainty on the strength of 
this limit. Black stars indicate the benchmark points chosen from each scenario, and they also indicate the region 

of parameter space that is not excluded by the other searches. 


ground. We see good agreement across all kinematic 
variables. In particular, unlike for other viable mod¬ 
els which tend to peak at high number of jets (large 
numbers of jets is often accompanied by large Ht), we 
can roughly reproduce the jet multiplicity plot distri¬ 
butions. If the excess persists this could be a powerful 


variable to discriminate between candidate interpreta¬ 
tions. We also note that this signal peaks at values for 
Ht and E™ lss far below the the thresholds for the kine¬ 
matic cuts of the ATLAS search, as opposed to mod¬ 
els based on the cascade decays of much heavier parti¬ 
cles that have been previously considered. Nonetheless, 










20 


> 

<L> 


O 


o 

o 


w 




Jet Multiplicity 


FIG. 4: The kinematic distributions compared to those of ATLAS in the on-Z search. The simulation predictions 
are summed with the ATLAS-calculated SM background to produce the model predictions. 


we still evade the bounds from the CMS search which 
has weaker E™ ISS cuts, due to the sharp increase in the 
background. For instance, the CMS background esti¬ 
mates were 478 ± 43, 39.2 ± 6.6, and 5.3 ± 2.3 events in 
the 100 GeV < Ef ss < 200 GeV, 200 GeV < £™ ss < 
300 GeV, and Ef iss > 300 GeV bins with n jets > 3, 
while for the f — C benchmark we predict 35, 11, and 
5.4 events respectively. We see that the large event 
rates in the low E™ ISS are well within the background 
uncertainties. 


B. Background contamination 


Another interesting feature of our signal is that it 
allows for the possibility of significant background con¬ 
tamination in the CMS search for the same final state 
described in section IIIB One of the most significant 


backgrounds in this search comes from SM Drell Yan 
(DY) production of Z bosons. To estimate this back¬ 
ground, the CMS collaboration used two independent 
data-driven methods and took a weighted average. One 
of these methods is based on the variable ‘jet -Z bal¬ 
ance’ (JZB) mm, which is important partic.ulatly in 
the high E™ lss search regions which constrain our sig¬ 


nal. The JZB of an event is defined by 

JZB = [ Pt\ - |^t Z) | = |i?T isS +Pt | - W ] \- 

jets 

(IV.l) 

SM processes like DY production typically result in JZB 
distributions that are symmetric about JZB = 0 GeV 
(because a non-zero value arises from jet energy resolu¬ 
tion effects), while some BSM processes can have JZB 
distributions that are strongly skewed towards positive 
values. This is expected when the Z is emitted back- 
to-back with an invisible particle, e.g. in a decay chain 
ending in \2 —► X^Z. For this reason, the JZB method 
estimates the DY background by assuming all events 
with JZB < 0 GeV are produced by DY, and extrapo¬ 
lating this to positive JZB values under the assumption 
that DY production is JZB-symmetric. 

It is clear therefore that signals with symmetric JZB 
distributions would contaminate this background esti¬ 
mate, reducing the sensitivity of the CMS search. In 
Fig. § we plot the JZB distributions after applying 
the cuts for the CMS nj e ts > 3, mid E™ lss bin for our 
three benchmark points, as well as for a gluino produc¬ 
tion model with the decay g —> qqx 2 —► qqx i% ( we have 
chosen the parameters rrig = 950 GeV, m x o = 50 GeV). 
We find that the JZB distribution is highly sensitive 
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FIG. 5: JZB probability distributions for the 
benchmark points when compared with a typical 
gluino production scenario (g —> qqx 2 —t qqXi %)■ The 
gluino scenario has the most positively skewed JZB 
distribution, while the JUS has almost symmetrical 
JZB. 


to the ti~t .2 mass splitting in this model. For small 
splitting, the Z tends to be very soft resulting in pos¬ 
itive JZB. For large splitting, the hard Z can result 
in symmetric or even negatively skewed JZB distribu¬ 
tions. The ^-S benchmark point has 41% of events with 
JZB < 0 GeV, comparing with only 8% in the gluino 
model. This also highlights the potential for the JZB 
distribution to be used as a discriminating variable be¬ 
tween new physics explanations of this excess should it 
persist in the next run of the LHC, due to its sensitivity 
to where the Z is emitted. 


V. CONCLUSION 

Motivated by the recent 3 a excess reported by the 
ATLAS collaboration in a Z + jets + I?™ lss channel we 
have studied if it can potentially be explained in the 
context of a natural supersymmetric spectrum involv¬ 
ing light stops. Strong constraints on such scenarios 
have led us to a compressed spectrum featuring two 
light and mixed stops and a light LSP. We identified 
three possible scenarios, characterized by flavor con¬ 
serving or flavor violating decays of the lightest stop, 
and the splitting between the two stop masses. We have 
shown that in all three scenarios it is possible to pro¬ 
duce the excess within 2 a, while in J^-C and Fq-C we 
can reproduce the excess within lcr of the ATLAS mea¬ 
surement. While the scenarios should be taken as ex¬ 
amples, it is clear that possible interpolations between 
them are capable of addressing the excess and would 


retain the same general features. Such features are a 
light stop with 225 GeV < < 325 GeV, almost de¬ 

generate with a Bino-like LSP and mixed with a second 
light stop with mass 325 GeV < < 550 GeV. 

The topology of the process differs from previous at¬ 
tempts to address the excess. The most substantial dif¬ 
ference is the production of the Z 's in the first step of 
a decay chain, and not in the last step in association 
with an invisible particle responsible for the E™ lss . In¬ 
terestingly, we have shown that it could lead to the con¬ 
tamination of background estimation based on the JZB 
method. This method is employed in a CMS search for 
a similar final state. We have estimated that as much 
as half of the signal could fall in the background control 
region, which could lead to over-exclusions. Addition¬ 
ally, we notice that the JZB variable could be used to 
discriminate between different signal topologies if this 
excess turns out to be due to new physics. 

Should this excess persist in run-II, it will be cru¬ 
cial to distinguish between the signal hypotheses. The 
signature proposed in this work is distinguished by its 
light compressed spectrum. This resulted in monojet 
searches being a highly sensitive probe of our signal. In 
addition, the search for b± —> 6y° is highly complemen¬ 
tary, and between these searches the region of parame¬ 
ter space which can explain the excess should be fully 
explored at 13 TeV. 

We note in passing that there are additional modest 
excesses of around two sigma or more in final states con¬ 
taining b-jets, leptons and MET, including a 1.9 sigma 
‘on -Z' excess in events with low jet multiplicity [48] , 
and various hints of same-sign dileptons with b-jets and 
MET (see [HJj for a summary). Light stops and sbot- 
toms can give rise to all of these signatures, and it is 
interesting to consider the possibility that if these really 
are all hints of new physics, they could have a unified 
explanation in a more complete model. Whether the 
ATLAS excess is a fluctuation or a first tantalizing hint 
of new physics will soon be decided. 
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